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1. Introduction 

Studying the statistical fluctuations of a system is a powerful method to characterise the ther- 
modynamic properties of a system. As a matter of fact, the presence of a phase transition is sig- 
nalled by an enhancement of the fluctuations in the system. The observables appropriate to extract 
this kind of information are called susceptibilities. Properly, the susceptibility describes the re- 
sponse of a system to an applied field. Given an operator we can usually determine a corresponding 
susceptibility that is the second derivative of the free energy density. In this paper we are interested 
in studying the quark number susceptibility (QNS) % = ^ : it is the response of n q , the quark 
number density, to an infinitesimal change in the quark chemical potential fx. 

The QNS has been studied so far at finite temperature and zero chemical potential, see Refs. [jj], 
|2[ |3j. This is largely because QNS is directly related to experimental measurements of fluctuations 
observed in heavy ion collisions [Eh. At non-zero density and zero (or low) temperature regime, 
QNS is not commonly studied; an application to QCD of rainbow approximation of the Dyson- 
Schwinger approach to study the QNS can be found in Ref. |5|]. One reason for the lack of studies 
in this regime is the well known sign problem: it is not possible, using standard tools, to simulate 
lattice QCD at non-zero density. 

In this paper we study the QNS at non-zero density and at low temperature in the context of 
SU (2) gauge theory, i.e. where simulations are feasible. 

2. Two color QCD 

In two color QCD, i.e. SU (2) gauge theory, quarks and antiquarks live in equivalent represen- 
tations of the group; the physical consequence is that qq mesons and qq, qq baryons are contained 
in the same hadron multiplet. In the limit when m K <C m p , i.e. when the pion mass is very small 
compared with the first non-Goldstone hadron, it is possible to study the system using the chi- 
ral perturbation theory (#PT) limit The fundamental result is that only for }X > jUo = \m n 
the quark number density n q becomes different from zero and at the same onset value /lo also a 
condensate (qq) ^ develops, signalling the spontaneous breakdown of the global £7(1) baryon 
number symmetry: a superfluid phase appears. In this phase there are tightly bound scalar diquarks 
but the hadrons are weakly interacting between them: therefore, just above the onset, the system is 
very dilute and is described as a Bose Einstein Condensate (BEC). 

Using #PT it is possible to determine the behaviour of different observables; in particular, we 
can write down the prediction for n q and its susceptibility %, at zero temperature T = and in the 
limit /I — > Ho (right-hand limit): 



where F is the pion decay constant and Nf the number of fermions. The diquark condensate, 
calculated in the same limits, is given in terms of the chiral condensate at zero chemical potential: 



n 



f T « 32N f F 2 (ji-fio), X xPT « ^N f F 2 - \92N f — (ji - juq) 



(2.1) 




(2.2) 
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There are reasons to think that above a certain value of the chemical potential, ji > \Lq, a de- 
generate system of weakly interacting quarks is more stable [0|. For an ideal gas (Stefan Boltzmann 
(SB) limit) of massless quarks and gluons, at T = 0, we have: 

N f N r 



™ NfN c , 



X 



SB 



(2.3) 



3n 2 • •- 

where N c is the number of colours. In this case the superfluidity is explained by a BCS condensation 
of Cooper pairs within a layer of thickness A around the Fermi surface; the diquark condensate is 
therefore given by: 

(qq) oc M 2 A . (2.4) 

Comparing Eq.s (|Q|), ( p72| ) with Eq.s ( |23| ), ( p~4| ), we see clearly that the two phases are charac- 
terised by two quite different behaviours. 

We consider also the order parameter related to the confinement property of the theory: the 
Polyakov loop L; as discussed in Ref. [||] the theory becomes deconfined only after }X > }Xd- Sur- 
prisingly, there is a regime, \Iq < ;li < /i£>, where the theory is confined, i.e. (L) = 0, but the other 
observables seem to suggest non interacting fermions: a confined BCS phase. This phase could be 
the so called quarkyonic phase introduced in Ref. [^] ; arguments against the extension of this idea 
to the case of N c = 2 can be found in Ref.[|l0|]. 



3. Calculation of the observables 



The fermion action with Nf = 2 and with a diquark source term J, necessary to study the 
diquark condensate, is given by: 



S f = yiM(h)yi + y 2 M{n)y 2 - 7y/i(Cy 5 )T 2 yf +^V2 r ( c 75)Wi 



(3.1) 



where M(jii) is the standard Wilson fermion matrix at non-zero chemical potential and C is the 
charge conjugation operator. If we introduce the change of variables = — y/ifCfy, <j> = C" X X%^2, 
Y = \ff\, Y = xjfi, it is possible to rewrite the action as: 



(#) 



-Jy 5 M(-ju) 




(3.2) 



It is worth mentioning that det {jjfl \M) = [det (M^M + Z 2 )] , therefore we can take the square 
root analytically, i.e. there is no square root problem. The partition function becomes: 



dUdVdVe- 5 *-*^ . 



It is now easy to write down an expression for n q using the matrix 



T dlnZ _ T 
V s ~diT~V s 



E<-*o 



d» J 



afi 



(Tr 



(3.3) 



(3.4) 



Moreover, from the definition of QNS we have: 



X 



djj. 



) + ( 



(3.5) 
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From this equation we can identify four different terms: 



71 



) 2 = -(Tr 
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T2 = +( 
CI = +( 
T3 = +( 
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(Tr 







Jt 
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(Tr 



(Tr 



-l 



>■ 



(3.6) 
(3.7) 
(3.8) 
(3.9) 



We see that from the second term of Eq. ( p.5| ) we get two terms, namely T2 and CI, because 
there are two ways to contract the spinors. 

The calculation of the traces is done by unbiased estimators, introducing complex noise 
vectors Y\ with the properties: (tj x ) = and (TfeTfy) = 8 xy . For example, the determination of the 
following trace, used for Tl and T2, is based on the relation: 



Tr 



1 



(3.10) 



Because for T2 we need two independent source vectors we refer to this term as disconnected term; 
the other three terms need only one source vector and we call them connected terms. 



4. Some numerical issues 

The source vector used to determine the traces can characterised by different noise distribu- 
tions. The standard method is based on the introduction of gaussian complex noise vector but it 
is possible also to use a Z2 complex noise vector (see Ref. [ |Tl| ] where the potential advantage are 
discussed). In the Z2 case, the complex noise vectors r\ takes one of the four values {±l,±i}, 
chosen independently with equal probability. 

In Table [j] we present an example where we used only three noise vectors, with the following 
parameters: j8 = 1.70, fc = 0.178, j = 0.04, pL = 0.25, 8 3 x 16. From this simple analysis, we can 





Tl 


T2 


T3 


CI 


X 


Gaussian 


2.25(8.30)E-06 


5.00(5.74)E-05 


0.4010(27) 


-0.3681(107) 


3.1(1.0) 


z 2 


6.3(10.0)E-06 


7.4(49.9)E-06 


0.3994(32) 


-0.3655(66) 


3.17(52) 



Table 1: Gaussian vs Z2 noise vectors. 

see that there is no any particular advantage in using a Z2 noise, so in the following we continue to 
use only gaussian noise distributions. 

We also tried to see what happen when we increase the number of noise vectors; in Table || 
we present an example with the same parameters used above. From this result it is evident how 
increasing the number of noise vectors has a strong effect mainly on observables with a small 
value; the effect is very limited if the observable is clearly different from zero. In Fig. [I] (Left) we 
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noise vectors 


Tl 


T2 


T3 


CI 


X 


3 


2.25(8.30)E-06 


5.00(5.74)E-05 


0.4010(27) 


-0.3681(107) 


3.1(1.0) 


300 


1.50(17)E-05 


1.53(12)E-05 


0.4013(14) 


-0.3612(20) 


4.04(1 l)E-02 



Table 2: The effect of different number of noise vectors. 
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Figure 1: (Left) The four terms of the QNS are plotted against the chemical potential. (Right) The ratio 
between Tl and T2 is plotted against the chemical potential. 



plot the four terms (note the sign of Tl) vs the chemical potential. The connected contribution CI 
gives clearly an important contribution either at low and high values of the chemical potential; note 
the changing of sign around a\i 0.66. The terms Tl and T2 are equal in magnitude but with 
opposite sign, therefore their contribution cancels everywhere, Fig. [I] (Right); in other words the 
variance of the quantity in Eq. ( 3.10 ) is equal to zero. We see very strong fluctuations, apparently 
not compatible with one, for /I < juo, then smoothed fluctuations comparable with one until Hd and 
after this a stable ratio equal to one. 



5. Numerical results 

In this contribution we are going to show some results obtained for /$ = 1.9 and K = 0.168. 
From previous works, see Ref. [Jsj] and the references therein, we know that for these parameters 
m n = 0.68(1), therefore ajUo ~ 0.34; moreover, there are signals of a BEC phase for a\i < 0.45. 

We obtained interesting results comparing the QNS with the other observables; here we want 
to stress an interesting effect we have observed studying the system at different temperatures. Note 
that the results we are presenting should eventually be extrapolated to j = 0.00. 

In Fig. U (Left) we plot the ratio #/jU 2 , for three different temperatures, versus the chemical 



potential. For an ideal gas of quarks this ratio would be a constant, see Eq. (2.3), and we see that a 
plateau is actually present for a\i < 0.55; after this value we can see a sharp increase of the QNS. 
Moreover, it is evident from this plot that the QNS is quite independent from the temperature, i.e. 
we do not see any drastic deviations in the behaviour of the three curves increasing \i. This is in 
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contrast with the Polyakov loop behaviour that shows deconfinement for three different values of 
the chemical potential, correspondingly at the three different temperatures, Fig. || (Right). 




It is instructive to compare our lattice numerical results with the equations corresponding to 
Eq. (2.3) but taking in account the finite volume and the lattice discretisation. In Ref. [{7|], see 
Eq. (26), the expression for the quark number density n SBL for free Wilson fermions on the lattice 
is presented. It is then easy to obtain % SBL . 

In Fig. H we plot the ratio between the measured QNS and x SBL for two values of the fermion 
mass: in the determination of % SBL we have to fix a value for the mass of the free fermion; unfor- 
tunately we do not know this value, therefore we consider massless fermions (note that this is the 
same limit used in Eq. (23)) and a value of the order of m K /2. In this case we observe a different 
behaviour for a\i < 0.45, reminiscent of the BEC phase, followed by a flat region compatible with 
ratio one, i.e. the system is behaving as free fermions, and then again we see an increase for higher 
values of \i. These plots again confirm the above scenario: we do not see any abrupt change for 
QNS, for any of the values of \i, where instead the Polyakov loop becomes different from zero. 

Note that the QNS is often taken as an alternative signal for deconfinement in lattice studies 
of the thermal QCD transition i.e. there is a strong connection between % and Polyakov Loop, see 
Refs. [ 12, 13|, [14]]. Our results shows that this relation is not replicated at low temperature and high 
density. 



6. Conclusion 



In this work we have studied the quark number susceptibility at non-zero density and low 
temperature in the case of two color QCD with two flavours. We have studied the contribution of 
the different terms, connected and disconnected, which contribute to it and finally we have shown 
its behavior at finite temperature. Surprisingly, we have seen that in this case there is no relation 
between the quark number susceptibility and the Polyakov loop, i.e. in this context it cannot be 
used as a signal for deconfinement as it is usually done at finite temperature and zero quark number 
density. This observation could suggest that the deconfinement transition is not characterised by 
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Figure 3: The ratio between the measured QNS and the ideal value for lattice free fermions for two values 
of the fermion mass: m = 0.00 and m — 0.34. The vertical dashed lines mark the position of jXo. 

a liberation of additional degrees of freedom; if this phenomenon is exclusive to two color QCD 
or related to the non small ratio m n jm p = 0.80(1) or, somehow, connected to the presence of 
a quarkyonic phase, is under study. Clearly, there is much to learn about deconfinement from 
studying a new physical environment. 
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